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Abstract The present paper provides a description of the 
oscillation code FILOU, its main features, type of applica- 
tions it can be used for, and some representative solutions. 
The code is actively involved in CoRoT/ESTA exercises (this 
volume) for the preparation for the proper interpretation of 
space data from the CoRoT mission. Although CoRoT/ESTA 
exercises have been limited to the oscillations computations 
for non-rotating models, the main characteristic of FILOU is, 
however, the computation of radial and non-radial oscilla- 
tion frequencies in presence of rotation. In particular, FILOU 
calculates (in a perturbative approach) adiabatic oscillation 
frequencies corrected for the effects of rotation (up to the 
second order in the rotation rate) including near degeneracy 
effects. Furthermore, FILOU works with either a uniform ro- 
tation or a radial differential rotation profile (shellular rota- 
tion), feature which makes the code singular in the field. 

Keywords methods: numerical • stars: evolution ■ stars: 
general • stars: interiors • stars: oscillations (including 
pulsations) • stars: rotation • (stars: variables:) delta Scuti 
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1 Introduction 



to correct the oscillation frequencies for the effects of rota- 
tion. In particular, the inclusion of these corrections (up to 
the second order including near degeneracy effects) to the 
oscill ation code and its numerical tests was part of my PhD. 
work ( ISuarezl 12002) (hereafter S02). In that work oscillation 
computations were extended to the case of models including 
a radial differential rotation profile Q = £2(r), i.e., a radial- 
dependent differential rotation (the so-called shellular rota- 
tion). This last characteristic is, by now, unique, and makes 
this code singular in the field. 

Although FILOU is currently optimised for the study of 
the pulsational behaviour of intermediate-mass classical pul- 
sators, namely 8 Scuti stars and y Doradus stars, the code is 
of universal use. It has been used, for instance, to model indi- 
vidua l 8 Scuti stars li ke the well-known Altair dSuarez et all 
l2005h . or 29 Cygnus dCasas et al.ll2006h. as wel l as to study 
8 Scuti stars in open clust ers dSuarez et alll2002HFox Machado et all 
l2006l : ISuarez et all 120071) : Moreover, it has serv ed to model 
high- amplitude 8 Scuti stars dPoretti et alll2005|) and to anal - 
yse the effec t of rotation on Petersen diagrams dSuarez et"ail 
2006a, 200 7[). Furtherm ore, it is worth highlighting the work 
bv lSuarez et all d2006bl) (from now on SGM06) which takes 
advantage of the code's main feature, i.e., the computation 
of adiabatic oscillation in presence of shellular rotation, to 
analyse the effect of such type of rotation on adiabatic os- 
cillation of moderately-fast rotating 8 Scuti stars. Concern- 



Numerous oscillation codes currently provide oscillation modesing J Doradus stars, FILOU has participated, additionally to 



for polytropes and ID models representative of different kind 
of pulsating stars. The history of the oscillation code FILOU 
is particularly associated with 8 Scuti stars. Originally de- 
veloped by F. Tran Min h and L. Leon at O bservatoire de 
Paris-Meudon (see lTran Miim&Leonl 119951) . the code has 
undergone several modifications and improvements in order 
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other modelling works, in one of the most recent and promis- 
ing asteroseismic tool for the modelling of su ch stars, the 
Freque ncy Ratio Method (FRM ), developed by iMova et al] 
d2005l) and lSuarez et all d2005h . 

From the point of view of the numerics, FILOU solves 
full sets of ODE (Ordinary Differential Equations) in a B VP 
(Boundary Value Problem), using a combined Galerkine - 
B-splines method which enhances the numerical precision 
with which the oscillation frequencies are calculated. Fur- 
thermore, as explained in the following sections, it is pos- 
sible to easily modify numerous numerical parameters in 
order to adjust the calculation optimally for the required 
model, which makes of FILOU a highly versatile code. 
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2 The adiabatic oscillations equations and boundary 
conditions 

FILOU is mainly based on th e oscillations equations a n d their 
pert urbations develop ed in iDziembowski & Goodd d 19921) 
and lSoufi etal.ldl998l) . In S02 we describe the second-order 
perturbation formalism used, which includes the effects of 
near degeneracy and considers the presence of a radial dif- 
ferential rotation (shellular rotation) profile, as well as its 
implementation in FILOU. 

The notation followed is similar to that used by many 
other oscillation codes, but adapted to the theoretical de- 
velopment considered. Although the code works with dif- 
ferent calculations schemes, namely, no rotation, Cowling 
approximation, and rotation (uniform and differential), in 
the present document only the most general case is con- 
sidered, i.e., the presence of shellular rotation. In such a 
case, oscillations are computed from the so-called pseudo- 
rotating models, which, as explained in S02, are constructed 
by modifying the stellar structure equations such as to in- 
clude the spherical symmetric contribution of the centrifu- 
gal acceleration, by means of an effective gravity g e g = g — 
A c {r) where g and A c (f) are the local gravity component 
and the centrifugal acceleration, respectively. The effects of 
the non-spherical components of the deformation of the star 
are included through a perturbation in the oscillation equa- 
tions. For instance, the perturbation of the mean density of 
a pseudo-rotating model po is considered of the form p 2 = 
P2i{r) P 2 (cos 0), where p 2 2(r) is denned in SGM06 (Eq. 15). 

Furthermore, when near degeneracy is taken into account, 
the eigenfrequency and the eigenfunction of a near-degenerate 
mode are then assumed of the form: 

co d = (Bo + 6)i + S>2 (1) 
k = I ajfaj + hj). (2) 

j=a,b 

where a>o = (dOo.a + ft)o.fo)/2 and a,- represent the coeffi- 
cients of the linear combination between the two con- 
sidered degenerate modes. Subscripts a and b represent 
whatever two rotationally coupled modes. First- and second- 
order corrections to the eigenfrequency in presence of near 
degeneracy are represented by 6)\ and ©2 respectively; £0./ 
and §i j are the non perturbed and first-order (see definitions 
and details in S02 and SGM06). The computation of individ- 
ual (Do./ as well as the corresponding zeroth- and first-order 
eigenfunctions is described in the next sections. 



where g e ff represents the effective gravity defined in the pre- 
vious section. The quantities pi and </>/, represent the eu- 
lerian perturbation of the pressure and the gravitational 
potential (definitions of individual terms can be found in 

FILOU and SGM06). Considering a differential rotation pro- 
file of the form: 



X2(r) = fl[l + ijo(r)] ! 



(5) 



where £2 represents the rotation frequency at the stellar sur- 
face, the eigenfrequencies (zeroth order) of a pseudo-rotating 
model are calculated from the linearised eigenvalue system: 
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which is solved by FILOU using the dimensionless variable 
x = r/R and R the stellar radius. In the calculation the fre- 
quency is expressed in terms of the dimensionless squared 
frequency 

,2 
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as well as other adiabatic quantities (see Appendix in SGM06) 
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As in SGD98, the following variables are also employed 
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rJ 
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(11) 



), Ux=U+ X . (12) 
A = Vg(joi ~yo2+yo3) A =£(£+!), (13) 



geff V dr 



where M and m r are the stellar mass and the mass enclosed 
in the sphere of radius r, respectively. 



2.1 Oscillation frequencies of a pseudo-rotating model 



In order to compute the oscillation frequencies of a pseudo- 
rotating model, (Oo.j the following dimensionless quantities 
are used 



yo\ = — , 
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yo4 
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2.2 The boundary conditions 

The system above (Eqs. [6]l is solved with the appropriate 
boundary conditions 



yo2 
yo\ - 



3yoi+}'04 = {1 = 0) 



^02^=0, y O 4-fy03=0 (t^O) 



"3 ; oiTr=0, 
j_ 



(14) 
(15) 
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Fig. 1 Normalised eigenfunctions yoi, yoi, yo3, and >>04, as a function of the normalised radial distance r/R, corresponding to the oscillation 
mode (n = S,£ = 1) calculated from a 1.8M@ 5 Scuti star model, with a rotational velocity of 100 kmr' at the stellar surface. 



at the centre of the star and, 
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(16) 
(17) 



(18) 



at the stellar surface. Figure [TJ illustrates the solutions for 
the normalised eigenfunctions _yoi, yo2, yo3, and yo4, corre- 
sponding to a non-radial mixed mode (n = 8, I = 1) obtained 
from a 1.8M Q model with a surface rotational velocity of 
lOOkms -1 , constructed with 2100 mesh points. 



2.3 First-order perturbed eigenfunctions 

When near-degeneracy effects are considered, first-order cor- 
rections to the eigenfunctions are required (see Eq.0. Con- 
sidering dimensionless variables equivalent to Eq.[3|-[4]with 
first-order perturbed quantities (£,\. r , </>{, p[), and the zeroth- 
order solutions obtained from Eq. [6] FILOU calculates such 



first-order perturbed eigenfunctions solving the following sys- 
tem: 



dyi , a A 
ax C r OQ 



(19) 



+ (yoi +Zo)(l + rj ) - (r]o + ai)Az 
= (C,a 2 - A*)yi + (A* + 1 - t^)ja - A*y 3 
<Tl + TJo)yol - (1 + T?0>0 



x— = (1 -U x )y 3 +y 4 



dx 

dn = 

fifoc 1 a,- 



(20) 



U 



A*y\ + V g y 2 - V g y3 + Ay 3 - U x y 4 



where X\ = V g (y i — y % +y$ ) and zo = y 02 /C^o • The horizon- 
tal component of £ j can be written as follows: 



yi , i + J]o 
zi = H : — yoi 



Col 



A 



1 + r?o 
A 



(21) 



where <7\ =Cl— Jo represents the first-order correction of 
the corresponding eigenfrequency. 
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2.4 The first- and second-order frequency corrections 

Second-order frequency corrections in presence of near de- 
generacy (Eq. Q]) are coded in FILOU using the following 
equations (S02, SGM06) 



in the case that S(Oo = COa, a — <Boi> 1S 0(Q ), and 



2,ab ' 
2\ 



(02 = 



8co? 



1\ 



8d)o / 



±\ "K, 



r(!) 



(22) 



(23) 



if 5(Oq is 0(12). In both cases first-order near degeneracy ef- 
fects are implicit, which are also calculated by FILOU using 



CO] = 



(Oih 



l.ab ; 



(24) 



Definitions of all terms involved are given in SGM06. The 
V and /I variables contains the corrections (up to the second 
order) for the effect of rotation on individual eigenfrequen- 
cies (Oq j obtained from Eq.[6] Such corrections are coded in 
FILOU using the Saio's notation, 



CO 



'7,2 



= eooj + n (C L - 1 -Jo) + — (d + m 2 D x \ 



(25) 



where Do = X\ +X2 and D\ = Y\ +Y2. Definitions and details 
can be found in (SGM06, Eqs 8-15). 



3 Structure & computation schemes 

FILOU is composed by a main program and some modules 
written in C, and two subroutines written in FORTRAN (77, 
95), which read input data from the equilibrium models and 
calculate the near degeneracy effects for the rotationally cou- 
pled modes. 

Computation of radial and non-radial oscillation frequen- 
cies of a given resonant cavity (input equilibrium model) is 
divided into three sequential steps: first, zeroth-order oscilla- 
tion frequencies (eigenvalue, (Oq) are computed as described 
in Sect. I2.lt then, for each eigenfrequency, the correspond- 
ing second-order frequency corrections (without including 
near degeneracy effects) are calculated (see Sect. |2.4| l. Fi- 
nally, the code selects, following certain rules (see S02 or 
SGM06), the rotationally coupled modes (only pairs of cou- 
pled modes are considered) and calculates their correspond 



L-Q 
1=1 
L=2 



Fig. 2 Effect of the number of mesh points on the adiabatic oscillation 
spectrum of a typical 5 Scuti star model (1.8M e ) as a function of 
the radial order n. Curves represent the frequency differences Av = 
|V2uoo — v 60o| obtained with 600 and 2000 mesh points. 



CESAM-type models: v3.*, 
GENEVE-type models 



v4* v5.* and 2k 



The input parameters, which are set by the user in a text 
file (ASCII), are read by the code when executed. The main 
parameters are: 

• The input equilibrium model file. 

• Type of computation. This option allows the user to force 
some kind of computing regime (for instance, Cowling 
approximation, no rotation, uniform rotation, differential 
rotation). As well, the user can choose the type of output 
files required, for instance (only the list of frequencies, 
include or not the corrections for the effect of rotation, 
near degeneracy effects, or even the eigenfunctions). 

• Frequency domain and spherical degree i range. 

• Type of boundary conditions (finite/infinite V, see Sect l2.ll) 

• Type of node assignation (zeros of yoi or JCD method). 

The basic outputs provided by FILOU are the list of eigen- 
frequencies and eigenfunctions. However, it is possible to 
obtain output files containing intermediate calculation data. 



4 Numerical techniques 



The numerical technique followed by FILOU is based on 
the Galerkin method, together with a finite sequence of B- 
splines, which is characterised for its flexibility, efficiency 
ing near degeneracy correcting terms (also described in Sect. l2.4jt id robustness. Although the code was conceived to solve 

the numerical problem of stellar non-radial oscillations, it 
actually provides solution to any non-linear system of func- 
tional equations, and covers several specific cases, such as 
the method of finite elements, Lagrangian and/or hermititian 
of any order, or even the Crank-Nicholson method of finite 
differences. 

Systems[6]and[20]are solved by approximating the eigen- 
functions (Eqs. [3]|4]l with B-Spline functions. The order of 



3.1 FILOU inputs & outputs 

FILOU inpunts are essentially some physical quantities (see 
Sect. 12.11) which are read from the equilibrium model and 
some initial parameters. Currently, the most updated version 
of the code allows the use of the following input models: 
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Fig. 3 Illustration of the impact of applying Richardson extrapolation 
on the oscillation frequency spectrum computations for the two models 
used in Fig. [2] The effect is shown through the frequency differences 
Av = |vre — v|, where Vre represents the oscillation frequencies ob- 
tained applying the Richardson extrapolation. 

such B-Splines functions can be chosen by the user, although 
optimum results are obtained typically for an optimal order 
between 4 and 6. The coefficients for each function are com- 
puted by integration following the Galerkin method. Firstly, 
the code scan the user-specified frequency range to obtain 
a first guess for the eigenfrequencies. Then, exact eigen- 
frequencies are searched using either the technique of di- 
chotomy or the Newton-Raphson method. 

It is worth highlighting the numerical versatility of the 
code, which can be optimised for the oscillations computa- 
tion of very different pulsating stars, g.e., solar-like pulsators 
(implying high-order p modes), g-mode pulsators (y Do- 
radus stars, white dwarfs), 8 Scuti stars, etc. This is so due to 
the numerous numerical parameters that can be adjusted. To 
name a few, the precision of the solutions (zeroth- and first- 
order eigenfrequencies and eigenfunctions) required, the size 
of the internal frequency interval in which the eigenfrequen- 
cies are searched for (this optimises the calculations in the 
cases of high-order and low-order frequencies), the thresh- 
old for valid solutions, etc. 



4. 1 Numerical tests & results 

In this section we report succinctly some of the numerical 
tests carried out on the code for its optimisation, as well as 
some of the tests carried out in the framework of the ESTA 
exercises. 

As has been shown during the ESTA exercises some of 
the most limiting aspects, from the numerical point of view, 
in the calculation of adiabatic oscillations are: the number 
of mesh points (including Richardson extrapolation), their 
distribution (rezonning), and the type of boundary condi- 
tions used. The current version of FILOU does not take nei- 
ther Richardson extrapolation nor rezonning into account 
(included in the next release of the code). Nevertheless we 
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Fig. 4 Illustration of the relative difference in frequency obtained, for a 
given main-sequence, 8 Scuti star model (1.8M ), w hen using V —> 
or V —> °° in the outer boundary conditions (see Sect. l2.lt . 

have carried out some numerical tests of such effects, which 
are illustrated in Figs. [2}{4l for a typical main-sequence, 5 
Scuti star model (~ 1.8 M Q ). Notice that all the effects cause 
a shift in the oscillation frequency which increases as far 
as the radial order increases. This behaviour is exponential 
when varying the number of mesh points (Fig. |2]l and when 
using Richardson extrapolation approximation (Fig. [3J. In 
both cases, for high-order p modes, the effects can reach up 
10 and 14;UHz, respectively. For the frequency domain of 8 
Scuti stars, i.e. n < 15, such effects can reach up 1 and IfiHz 
which are nonnegligible compared with the high-precision 
in frequency detection that CoRoT mission is expected to 
provide. 



4.2 Results. The effect of shellular rotation on adiabatic 
oscillations 

One of the most important works carried out using the oscil- 
lation code FILOU is the study of the effect of shellular ro- 
tation on adiabatic oscillations (SGM06). Indeed, that work 
takes advantage of the main feature of the code, i.e., the cal- 
culation of oscillations in presence of a radial differential ro- 
tation. Figure |5] (left panel) illustrates the effect of a shellular 
rotation on the radial displacement eigenfunction you for a 
mixed mode obtained for a 1 .8M , 8 Scuti star model. Such 
eigenfunctions can be obtained with FILOU when calculating 
the oscillation spectra of pseudo-rotating models (see S02) 
computed assuming local conservation of the angular mo- 
mentum. Furthermore, the effect of shellular rotation on the 
oscillation frequencies is also significant for 8 Scuti stars. 
In Fig. HJ such an effect is depicted as a function of the 
radial order. In that figure A a> d and Aco represent mode- 
to-mode frequency differences with (bottom panel) and 
without (top panel) taking the effect of near degener- 
acy into account. These quantities are calculated, for a 
given mode, as the difference between the oscillation fre- 
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mixed mode n=2, 1=1 
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Fig. 5 Left panel displays weighted radial displacement eigenfunctions for a mixed mode as a function of the normalised radial distance r/R. 
Solid and dashed lines represent the / function computed for a shellular-rotating model. Dash-dotted lines represent the / function computed for 
a uniformly-rotating model. Dotted lines represent the rotation profile (scaled in the figure for clarity) given by the radial function T}o(r). Right 
panel shows mode-to-mode frequency differences between differentially and uniformly-rotating 1.8M© models. Symmetric, solid-line branches 
represent from top to bottom, differences for m = — 1 and m = + 1 mode frequencies respectively. For m = modes, differences are represented 
by a dotted line. The shaded region represents an indicative frontier between the region of g and gp modes (left side) and p modes (right side). 
Taken from SGM06. 



quency obtained assuming shellular rotation with the os- 
cillation frequency obtained assuming uniform rotation. 
The term Qj represents the additional effect of near de- 
generacy on the oscillation frequencies for two degener- 
ate modes and j (see SGM06 for more details). As shown 
in SGM06, for g and mixed modes, significant effects (up 
to 3 juHz) on the the oscillation frequencies are predicted. 
For high-frequency p modes, such effects can reach up to 
1 jUHz. Such effects are likely to be detectable with CoRoT 
data, provided numerical eigenfrequencies reach the level of 
precision required. See SGM06 for a detailed discussion on 
these results. 
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